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Abstract 

Using data accumulated with the CLEO detector corresponding to an integrated luminosity of 
£=5.63 pb- 1 on the peak of the ip(2S) (3.08 x 10 6 ip(2S) decays) and 20.70 pb" 1 at ^/i=3.67 GeV, 
we report first measurements of the branching fractions for the following 13 decay modes of 
the ip(2S): 7?3vr, r/'3vr, P K + R-, K+R-tt+tt-tt , 2(K + K~), 2(K + K-)tt°, ppp, ppir+ir-ir , npp, 
ppK + K~, AA"7r + 7r~, ApK + , and ApK + ir + ir~, and more precise measurements of 8 previously 
measured modes: 2(ir + TT~), pir + ir~, 2(tt + tt^)tt° , ljit + it~, K + K~it + it~ , loK + K~, (J)K + K~, and 
ppTT + TT~ . We also report new branching fraction measurements of c/>7r + 7r~ and ujpp and upper limits 
for r/7r + 7r~, i]K + K~ and <ppp. Results are compared, where possible, with the corresponding J/ijj 
branching ratios to provide new tests of the 12% rule. 
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The states J/if> and ip{2S) are non-relativistic bound states of a charm and an anti-charm 
quark. In perturbative QCD the decays of these states are expected to be dominated by the 
annihilation of the constituent cc into three gluons or a virtual photon. The partial width 
for the decays into an exclusive hadronic state h is expected to be proportional to the square 
of the cc wave function overlap at zero quark separation, which is well determined from the 
leptonic width 1.]. Since the strong coupling constant, a s , is not very different at the J /if) 
and if)(2S) masses, it is expected that for any state h the J /if) and if)(2S) branching ratios 
are related by |2| 

_ sms) -> h) „ isms) -> en-) _ (l27±(] , )% m 

where B denotes a branching fraction, and the leptonic branching fractions are taken from the 
Particle Data Group (PDG) This relation is sometimes called "the 12% rule". Modest 
deviations from the rule are expected [3] • Although the rule works well for some specific decay 
modes of the if>(2S), isospin conserving ip(2S) decays to two-body final states consisting of 
one vector and one pseudsoscalar meson exhibit strong suppression: jm is suppressed by 
a factor of seventy compared to the expectations of the rule [H 0, IE El- Also, vector- 
tensor channels such as pa2(1320), and K* (892)^2 (1430) are significantly suppressed [H, @|. 
A recent review of relevant theory and experiment concludes that current theoretical 
explanations are unsatisfactory. Clearly, more experimental results are desirable. 

This Letter presents new measurements of a wide selection of if)(2S) decays including 
modes with and without strange particles and with and without baryons. The following 
modes of the ip{2S) are observed for the first time: t]3tt, i]'3tt, pK + K~, K + K~ tt + tt° , 
2(K + K~), 2(K + K-)tt°, ppp, ppn + n-ir° , rjpp, ppK + K~, AA7r + 7r~, ApK + , and ApK+n+n-. 
More precise branching fraction measurements of previously measured modes 2(7r + 7r~), 
p7r + 7r _ , 2(7T + 7r~)7r°, c<J7r + 7r~, K + K--k + -k~ , u)K + K~, <f>K + K~, and ppiT + 7f- are also re- 
ported. We also measure 07r + 7r~ and uopp and obtain upper limits for i]tt + tt~, t]K + K- 
and (ppp. Where applicable, the inclusion of charge conjugate modes is implied. As 14 of 
the modes we study have been previously observed at the J /if), 14 tests of the 12% rule are 
made, 5 for the first time and 9 with greater precision than corresponding previous tests. 

The data sample used in this analysis is obtained at the ip(2S) and the nearby continuum 
in e + e~ collisions produced by the Cornell Electron Storage Ring (CESR) and acquired 
with the CLEO detector. The CLEO III detector features a solid angle coverage for 
charged and neutral particles of 93%. The charged particle tracking system, operating in 
a 1.0 T magnetic field along the beam axis, achieves a momentum resolution of ~0.6% at 
p — 1 GeV/ c. The calorimeter attains a photon energy resolution of 2.2% at E 1 — 1 GeV and 
5% at 100 MeV. Two particle identification systems, one based on energy loss (dE/dx) in the 
drift chamber and the other a ring-imaging Cherenkov (RICH) detector, are used together 
to separate K ± from n^. The combined dE/dx-RlCH. particle identification procedure has 
efficiencies exceeding 90% and misidentification rates below 5% for both tt^ 1 and K ± for the 
momenta below 2.5 GeV/c. 

Half of the if)(2S) data and all the continuum data were taken after a transition to the 
CLEO-c detector configuration, in which CLEO Ill's silicon-strip vertex detector was 
replaced with a six-layer all-stereo drift chamber. The two detector configurations also 
correspond to different accelerator lattices: the former with a single wiggler magnet and a 
center-of-mass energy spread of 1.5 MeV, the latter (CESR-c [9]) with six wiggler magnets 
and an energy spread of 2.3 MeV. 
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The integrated luminosity (£) of the datasets was measured using e + e~, 77, and 
final states [l0|. Event counts were normalized with a Monte Carlo (MC) simulation based 
on the Babayaga event generator combined with GEANT-based jl^ detector modeling. 
The data consist of £=5.63 pb" 1 on the peak of the ip(2S) (2.74 pb" 1 for CLEO III, 2.89 pb" 1 
for CLEO-c) and 20.70 pb" 1 at ^=3.67 GeV (16 MeV below the i/>{2S), all CLEO-c). 
The nominal scale factor used to normalize continuum yields to ip(2S) data is f nom = 
0.270 ± 0.005, and is determined from the integrated luminosities of the data sets corrected 
for the 1/s dependence of the cross section, where the error is from the relative luminosity 
uncertainty, and the uncertainty in the s dependence of the cross section. (The scale factor 
does not differ appreciably if a high power of 1/s is used: a factor of 0.87% for each power). 
The value of / used for each mode also corrects for the small difference in efficiency between 
the ip(2S) and continuum data. 

Standard requirements are used to select charged particles reconstructed in the tracking 
system and photon candidates in the Csl calorimeter. We require tracks of charged particles 
to have momenta p > 100 MeV/c and to satisfy | cos#| < 0.90, where 9 is the polar angle 
with respect to the e + direction. Each photon candidate satisfies E 1 > 30 MeV and is more 
than 8 cm away from the projections of tracks into the calorimeter. Particle identification 
from dE/dx and the RICH detector is used on all charged particle candidates. Pions, kaons, 
and protons must be positively and uniquely identified. That is: pion candidates must 
not satisfy kaon or proton selection criteria, and kaon and proton candidates obey similar 
requirements. Charged particles must not be identified as electrons using criteria based on 
momentum, calorimeter energy deposition, and dE/dx. 

The invariant mass of the decay products from the following particles must lie within lim- 
its determined from MC studies: vr° (120 < M 77 < 150 MeV), 77 (500 < M 77 < 580 MeV), 
77 (530 < M n+7r - w o < 565 MeV), u (740 < M w+rf < 820 MeV [760 < M T+n - n o < 800 MeV 
for the upp final state]), <p (1.00 < M K+K - < 1.04 GeV), and A (1.112 < M pn - < 
1.120 GeV). For 7r° — > 77 and 77 — > 77 candidates in events with more than two pho- 
tons, the combination giving a mass closest to the known 7T° or 77 mass is chosen, and a 
kinematically constrained fit to the known parent mass is made. Fake 7r°'s and rj's are sup- 
pressed by requiring that each electromagnetic shower profile be consistent with that of a 
photon. For 77 — ► 7r + 7r~7r° and uj — > 7r + 7r~7r°, the 7r° is selected as described above, and 
then combined with all possible combinations of two oppositely charged pions choosing the 
combination that is closest to the rj(u) mass. A kinematically constrained fit is not used for 
either of these modes, or for <f/s or for A's. For A — > pir~ , a fit of the pn~ trajectories to a 
common vertex separated from the e + e~ interaction ellipsoid is made. Contamination from 
Kg decays is eliminated by the energy and momentum requirements imposed on the event, 
and by particle identification. 

Energy and momentum conservation requirements are imposed on the reconstructed final 
state hadrons, which have momentum pi and combined measured energy E v - ls . We require the 
measured scaled energy E vis /E c ^ be consistent with unity within experimental resolution, 
which varies by final state. We require | J2Pi\/E c .m. < 0.02. Together these requirements 
suppress backgrounds with missing energy or incorrect mass assignments. The experimental 
resolutions are smaller than 1% in scaled energy and 2% in scaled momentum. 

For the final states with four charged tracks and a 7r°, an additional cut is applied to 
remove a background of radiative events. When the photon is combined with a low-energy 
photon candidate, it can imitate a tt°. We require (-EWacks + Ery)/E c _ m _ < 0.995, where E 1 
is the energy of the highest energy photon. 
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In order to compute Qh in modes with two or more charged pions, two 7r 0/ s or an T], it is 
necessary to remove the contribution from the transitions ip(2S) — > J/ipX, where X = 7r + 7r~, 
7r°7r°, or rj ^3[. Accordingly we reject events in which the mass of any of the following falls 
within the range 3.05 < m < 3.15 GeV 0]: the two highest momentum oppositely charged 
tracks, the recoil mass against the two lowest momentum oppositely charged tracks, or the 
mass recoiling against the 27r°'s or 77. 

For every final state, a signal selection range in E Y i B / E c _ m _ is determined by Monte Carlo 
simulation, and a sideband selection range is defined to measure background. Final states 
with an intermediate p, rj, rj', u, <fi or A must satisfy a scaled energy signal selection range 
requirement identical to the corresponding mode without the intermediate particle. For 
example, the scaled energy signal selection range is the same for <pK + K~ and K + K~ K + K~ . 
For final states with an rj, rj', u>, <p, or A the event yield is determined from signal and sideband 
selection ranges of the particle mass. For final states with a p the event yield is determined 
from a fit to the n + 7r~ invariant mass by a Breit-Wigner with parameters taken from the 
PDG i (for p7r + 7r , the same-sign 7r ± 7r ± invariant mass distribution is used to estimate 
the background); an exception is ppp, where the yield is determined from the signal region: 
0.54 < m n7T < 1.0 GeV assuming a linear background. Most modes studied in this Letter 
have resonant sub-modes, however; we only report branching ratios for resonant sub-modes 
which can be cleanly separated from background. 

In Fig. ^ the scaled energy and invariant mass distributions are shown for three typical 
modes: 2(tt + tt~), 2(tt + 7t~)tt° , and K + K~tt + tt~ . Evidence for production of the p, u and 4> 
resonances, respectively, is observed in the corresponding mass spectra. 

Event totals are shown for both the ip(2S) and the continuum in Table EJ where S^ps) 
(S co ) is the number of events in the signal region and B^s) (B co ) the number of events 
in the sideband region in ip(2S) (continuum) data. Under the assumption that interference 
between ip(2S) decay and continuum production of the same final state is absent, the number 
of events attributable to each ip(2S) decay mode, Ns, is 

Ns = S^s) - B^2S) - f(S co - B co ), (2) 

where / is mode dependent and listed in Table |TJ We observe a statistically significant 
signal in all modes except 7]7t + tt~, r]K + K~ and (ppp. The signal with the least statistical 
significance is r\pp (4.2a). 

The efficiency, e, for each final state is the average obtained from MC simulations [12j for 
both detector configurations; the two values are typically within a few percent (relative) to 
each other. No initial state radiation is included in the Monte Carlo, but final state radiation 
is accounted for. The efficiencies in Table E include the branching ratios for intermediate 
final states. 

We correct Ns by the efficiency and normalize to the number of ip(2S) decays in the data: 
3.08 x 10 6 determined by the method described in The resulting branching ratios are 
listed in Table HI along with a comparison to the PDG pj. With the exception of 2(7r + 7r~) 
and p7T + 7T _ , none of the branching ratios in |l| were corrected for the contribution from 
continuum production. 

The systematic errors on the ratio of branching fractions share common contributions 
from the number of ip(2S) decays (3.0%), uncertainty in / (2.0%), trigger efficiency (1%), 
and Monte Carlo statistics (2%). Other sources vary by channel. We include the fol- 
lowing contributions for detector performance modeling quality: charged particle track- 
ing (0.7% per track), ■n Q / r q(-^ 77) finding (4.4%), A finding (3%), n/K/p identification 
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FIG. 1: Distributions for the modes 2(7t + 7t~) (a) and (b), 2(7r + 7r~)7r° (c) and (d), and 
K + K~it + it~ (e) and (f). The pairs of arrows indicate the signal regions, (a), (c), and (e) The 
scaled total energy (filled circle with error bar: ip(2S) data, solid line: Monte Carlo, shaded his- 
togram: continuum), (b) The tt + tt~ invariant mass in 2(tt + tt~) (solid line: data, dashed line: like 
sign it^ir^ invariant mass, which is used to estimate combinatorial background), (d) The tt + tt^tt° 
invariant mass in 2(7r + 7r~)7r (filled circle with error bar: ip(2S) data, solid line: Monte Carlo, 
shaded histogram: continuum), (f) The K + K~ invariant mass in K + K~n + ir~ . (The symbol key 
is the same as (d).) 



(0.3%/1.3%/1.3% per identified n/K/p), and scaled energy and mass resolutions (2%). The 
systematic error associated with the uncertainty in the level of background is obtained by 
recomputing the branching ratio when the background at the ip(2S) and the continuum 
are coherently increased by la. Since the background in many modes is small, the Poisson 
probability for the observed number of background events to fluctuate up to the 68% C.L. 
value is calculated and interpreted as the uncertainty in the level of background. Many of 
the modes studied have resonant submodes, for example: 2(7r + 7r~) where pnir is dominant, 
UJ717T where &i7r is dominant and 07/2(1270) is significant, K + K~7t + 7t~ where K*(892)Kn and 
pK + K~ are large and 07r + 7r~ is small, and r/37r where r)7r°p°, i]p + ir~ , and r/p~7r + dominate. 
The efficiencies for the modes with resonant submodes were weighted using the fractions of 
the submodes. Allowing for the presence of resonant submodes changes the efficiency by 
less than 5% relative to the non-resonant efficiency for most modes we have studied. The 
difference between the weighted efficiency and phase space efficiency is taken as a measure 
of the uncertainty in the efficiency. For those modes with no observed sub-modes, a conser- 
vative 10% uncertainty is assigned for decay model dependence. Systematic uncertainties 
are significant for all modes and the dominant error for many. The measurements in this 
Letter are in reasonable agreement with previous measurements where such exist p|. 

For the 14 modes where the same final state has been previously measured at the J/i/j, 
the Qh value is computed using the absolute ip(2S) branching ratios determined in this 
analysis and J / ip branching ratios from and given in Table U For 5 of these modes: 
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TABLE I: For each final state h the following quantities are given: the number of events in the 
signal region, S co , and background from sidebands, B co , in continuum data; the scale factor, /; the 
number of events in the signal region, Sw25), an d background from sidebands, B^s), in ipftS) 
data; the number of events attributable to ip(2S) decay, Ns, computed according to Eq. [21 the 
average efficiency, e; the absolute branching fraction with statistical (68% C.L.) and systematic 
errors; previous branching fraction measurements from the PDG [lj], an d the Qh value. For r/37r, 
the two decays modes r]3ir(r] — ► 77) and 77371- (77 — > 3ir) are combined on line r]3ir. 



mode 


Sco 


B co 


/ 




B 1 p(2S) 


N s 


£ 


B(ip(2S) -» h) 


B (PDG) 


Qh 


h 
















(units of 10" 4 ) 


(units of 10" 4 ) 


(%) 


2(ir+ir-) 


1471 


28 


0.2668 


713 


20 


308.0 


0.4507 


2.2±0.2±0.2 


4.50±1.00 


5.55±1.53 




1168 


- 


0.2667 


597 


- 


285.5 


0.4679 


2.0±0.2±0.4 


4.20±1.50 


- 


2(7r + vr-)7r 


352 


25 


0.2550 


1825 


39 


1702.6 


0.2115 


26.1±0.7±3.0 


30.00±8.00 


7.76±1.10 


r]ir + ir~ 


15 





0.2501 


13 


2 


7.2 


0.0416 


< 1.6 


- 


- 


LUir + ir~' 


43 


9 


0.2357 


437 


38 


391.0 


0.1553 


8.2±0.5±0.7 


4.80±0.90 


11.35±1.94 


r]3ir(ri — > 77) 


27 


2 


0.2513 


243 


35 


201.7 


0.0639 


10.3±0.8±1.4 


- 


- 


7]3ir(rj — > 37r) 


20 


9 


0.1820 


53 


1 


50.0 


0.0199 


8.1±1.4±1.6 


- 


- 


7?37T 
















9.5±0.7±1.5 


- 


- 


?/37r 


1 





0.1721 


17 


4 


12.8 


0.0092 


4.5±1.6±1.3 


- 


- 


ET+A'-Tr+Tr- 


871 


83 


0.2688 


1072 


43 


817.2 


0.3742 


7.1±0.3±0.4 


16.00±4.00 


9.85±3.23 


pK+K- 


170 


- 


0.2602 


268 


- 


223.8 


0.3361 


2.2±0.2±0.4 


- 


- 


(/>7T + 7r~ 


33 


13 


0.2703 


73 


20 


47.6 


0.1744 


0.9±0.2±0.1 


1.50±0.28 


11.07±3.30 


A+A-vr+vr-vr 


634 


18 


0.2556 


888 


19 


711.6 


0.1818 


12.7±0.5±1.0 


- 


10.59±2.81 


r]K + K- 


3 





0.2396 


7 


2 


4.3 


0.0354 


< 1.3 


- 


- 


ojK + K- 


62 


12 


0.2435 


97 


8 


76.8 


0.1288 


1.9±0.3±0.3 


1.50±0.40 


10.19±2.96 


2(K + K~) 


100 


11 


0.2669 


85 


2 


59.2 


0.3118 


0.6±0.1±0.1 




6.71±2.74 


0A+A- 


46 


15 


0.2642 


49 


4 


36.8 


0.1511 


0.8±0.2±0.1 


0.60±0.22 


5.14±1.53 


2(A + A-)vr 


20 





0.2675 


51 


1 


44.7 


0.1339 


1.1±0.2±0.2 








337 


28 


0.2509 


1010 


28 


904.5 


0.4943 


5.9±0.2±0.4 


8.00±2.00 


9.90±1.16 


ppp 


23 




0.2570 


67 




61.1 


0.4119 


0.5±0.1±0.2 








204 


9 


0.2312 


499 


19 


434.9 


0.1921 


7.3±0.4±0.6 




18.70±5.80 


7]pp 


2 


1 


0.2350 


12 


2 


9.8 


0.0399 


0.8±0.3±0.3 




3.80±2.09 


uopp 


26 


4 


0.2173 


37 


11 


21.2 


0.1129 


0.6±0.2±0.2 


0.80±0.32 


4.69±2.22 


ppK + K~ 


25 


1 


0.2478 


37 


1 


30.1 


0.3671 


0.3±0.1±0.0 






4>pp 


2 


3 


0.2631 


6 


2 


4.3 


0.1732 


< 0.24 


< 0.26 




AA7r + 7r~ 


23 


4 


0.1902 


91 


14 


73.4 


0.0844 


2.8±0.4±0.5 






ApA+ 


65 


7 


0.2586 


97 


8 


74.0 


0.2472 


1.0±0.1±0.1 




10.92±2.93 


ApA + 7r + 7r _ 


29 


3 


0.1631 


57 


7 


45.8 


0.0847 


1.8±0.3±0.3 







K + K^n + 'n-^7T , 2(K + K~), pp7r + 7r~7r°, rjpp and ApK + , this is the first measurement of 
Q h . All modes except pp7i + n^n° have Qh below 12.7%. The modes <fiK + K~, 2(n + n^), 
2(7r + 7r~)7r°, rjpp, and upp are suppressed with a significance of 4.7, 4.4, 4.1, 4.1, and 3.5 
combined statistical and systematic standard deviations, respectively. 
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In summary we have presented first branching fraction measurements for 13 modes of the 
i/;(2S): T]3tt, T]'3tt, pK + K~, K + K-n + n-iT , 2(K + K~), 2(K + K-)n°, ppp, ppvr+vr-vr , ripp, 
ppK + K~ , AAtt + tt~, ApK + , ApK + iT + 7T~ , and more precise measurements of 8 previously 
measured modes: 2(7r + 7r~), /wr + 7r _ , 2(7r + 7r~)7r°, ujtt + 7t~, K + K~tt + tt~ , uK + K~, <pK + K~, 
ppTT + Ti~ . We also measure the branching ratios for 07r + 7r~ and ujpp and obtain upper limits 
for i]7r + 7T~ , r]K + K~ and (ppp. A full set of measurements of tt + 7t~+X, K + K~ +X and pp+X 
where X is p/rj/u/cj) is included in this analysis. Results are compared, where possible, with 
the corresponding J/ip branching ratios to test the 12% rule. Five modes are inconsistent 
with the rule, but all values of Qh are within a factor of 2 of 12.7%. Since this analysis 
covers a wide variety of final states including modes with and without strange particles and 
with and without baryons, the pattern of branching ratios and Qh values provide tests of 
theoretical models and should allow new insight into the production widths and the final 
state interactions operative in the decays of the ip(2S). 
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